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Abstract 

We present a flavor model of quarks and leptons with the non-Abelian discrete symmetry 
54 in the framework of the SU{5) SUSY GUT. Three generations of 5-plets in SU{5) 
are assigned to 3 of 84 while the first and second generations of 10-plets in SU{5) are 
assigned to 2 of ^4, and the third generation of 10-plet is assigned to 1 of 5*4 . Right- 
handed neutrinos are also assigned to 2 for the first and second generations and 1' for the 
third generation. We predict the Cabibbo angle as well as the tri-bimaximal mixing of 
neutrino flavors. We also predict the non-vanishing Ues of the neutrino flavor mixing due 
to higher dimensional mass operators. Our predicted CKM mixing angles and the CP 
violation are consistent with experimental values. We also study SUSY breaking terms in 
the slepton sector. Our model leads to smaller values of flavor changing neutral currents 
than the present experimental bounds. 
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1 Introduction 



There are many free parameters in the standard model including its extension with neutrino 
mass terms and most of them are originated from the flavor sector, i.e. Yukawa couplings of 
quarks and leptons. Quark masses and mixing angles have been discussed in the standpoint 
of the flavor symmetries. The discovery of neutrino masses and the neutrino flavor mixing has 
stimulated the work of the flavor symmetries. Recent experiments of the neutrino oscillation 
go into a new phase of precise determination of mixing angles and mass squared differences 
[H El El S] , which indicate the tri-bimaximal mixing for three flavors in the lepton sector [3 
El m [S]- These large mixing angles are completely different from the quark mixing ones. 
Therefore, it is very important to find a natural model that leads to these mixing patterns 
of quarks and leptons with good accuracy. 

The flavor symmetry is expected to explain the mass spectrum and the mixing matrix of 
both quarks and leptons. Especially, the non-Abelian discrete symmetry [9] has been studied 
intensively in the quark and lepton sectors. Actually, the tri-bimaximal mixing of leptons has 
been at first understood based on the non-Abelian finite group ^4 [lOl [TTl [121 ESI E] • Until 
now, much progress has been made in the theoretical and phenomenological analysis of A4 
flavor model [T5]-[73]. On the other hand, much attention has been devoted to the question 
whether these models can be extended to describe the observed pattern of quark masses 
and mixing angles, and whether these can be made compatible with the SU{5) or 5*0(10) 
grand unified theory (GUT). The attractive candidate is the 6*4 symmetry, which has been 
already used for the neutrino masses and the neutrino flavor mixing [TH [751 [ZS [77]. The 
exact tri-bimaximal neutrino mixing is realized in 6*4 flavor models [781 [HOl [SI [821 [Ml [H] . 
Many detail studies in the 5*4 flavor model have been presented for the quark and lepton 
sectors [85]-|95]. There are attempts to unify the quark and lepton sectors toward a grand 
unified theory of flavor [86l [13, [HE] , however, quark mixing angles were not predicted clearly. 

Recently, 5*4 flavor models to unify quarks and leptons have been presented in the frame- 
work of the SU{5) SUSY GUT [80] or ^O(IO) SUSY GUT [96]. However, quantitative 
analyses have not been given there since the contribution from higher dimensional mass op- 
erators are not discussed. There also appeared the flavor model in SU{5) SUSY GUT [HZ] 
and the Pati-Salam SUSY GUT [HHI, taking account of higher dimensional mass operators. 

In this paper, we present another S4 flavor model with Z4 taking account of higher dimen- 
sional mass operators. We predict the deviation from the tri-bimaximal mixing of the lepton 
flavor numerically. The CKM mixing angles and CP violation are discussed numerically 
owing to higher dimensional mass operators. We also discuss the flavor changing neutral 
current (FCNC) in the SUSY sector, which is important to constrain the parameter of the 
flavor model. 

The 5*4 group has 24 distinct elements and irreducible representations 1, 1', 2, 3, and 
3'. Three generations of 5-plets in SU{5) are assigned to 3 of S4 while the first and second 
generations of 10-plets in SU{5) are assigned to 2 of 5*4, and the third generation of 10-plet 
is assigned to 1 of S4. These assignments of 6*4 for 5 and 10 lead to the completely different 
structure of quark and lepton mass matrices. Right-handed neutrinos, which are SU{5) 
gauge singlets, are also assigned to 2 for the first and second generations, and 1' for the third 
generation. These assignments realize the tri-bimaximal mixing of neutrino flavors. Gauge 
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singlet scalars, which are so called flavon, are introduced, in our model. Relevant vacuum 
alignment of flavons gives the quark flavor mixing angles as well as the tri-bimaximal mixing 
of neutrino flavors. Especially, the Cabibbo angle is predicted to be around 15°. 

In section 2, we present the 5*4 x Z4 x f/(l)^7v flavor model of quarks and leptons in SU{5) 
SUSY GUT, and discuss the effect of the higher dimensional mass operators. The deviation 
from the tri-bimaximal mixing is predicted, and CKM mixing angles and the CP violation are 
discussed in detail. In section 3, the alignment of the vacuum expectation values (VEVs) is 
derived. In section 4, the FCNC in the slepton sector is discussed. Section 5 is devoted to the 
summary. We present the multiplication rule of 5*4, the determination of U{1)fn quantum 
numbers, and the analysis of the scalar potential in Appendices A, B, and C, respectively. 

2 ^4 X Z4 X U{1)fn flavor model with SU{b) SUSY GUT 
2.1 Assignments of superfields 

We present the 5*4 flavor model in the framework of SU (5) SUSY GUT. The flavor symmetry 
of quarks and leptons is the discrete group 5*4 in our model. The group 5*4 has irreducible 
representations 1, 1', 2, 3, and 3'. The multiplication rule is shown in Appendix A. 
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Table 1: Assignments of SU{5), S4, Z4, and U{1)fn representations. 



Let us present the model of the quark and lepton flavor with SU (5) SUSY GUT. In SU (5), 
matter flelds are unifled into 10 and 5 dimensional representations. Three generations of 5, 
which are denoted by Fj (i = 1, 2, 3), are assigned to 3 of S4. On the other hand, the third 
generation of the 10-dimensional representation is assigned to 1 of S'4, so that the top quark 
Yukawa coupling is allowed in the tree level. While, the flrst and second generations are 
assigned to 2 of S'4. These 10-dimensional representations are denoted by T3 and (Ti,T2), 
respectively. Right-handed neutrinos, which are SU{5) gauge singlets, are also assigned to 
1' and 2 for A''^ and (A'g, A"^), respectively 0. 

We introduce new scalars Xi addition to the 5-dimensional, 5-dimensional, and 45- 
dimensional Higgs of SU{5), H^, ifg, and H4^, which are assigned to 1 of S'4. These new 

^ Our S'4 assignments of matter fields are same as ones in the model |97j except that right-handed neutrinos 
are assigned to 3 there. 
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scalars are supposed to be SU{5) gauge singlets. Scalars (X17X2) and (x3)X4) are assigned 
to 2, (X5,X6,X7) are assigned to 3', (x8,X9,Xio) and (xn, X12, Xis) are assigned to 3, and 
Xi4 is assigned to 1 of 5*4 representations, respectively. In the leading order, (xs^X-i) are 
coupled with the right-handed Majorana neutrino sector, (xs^Xe^X?) are coupled with the 
Dirac neutrino sector, (xs^Xq^Xio) and (xii; X12, X13) are coupled with the charged lepton 
and down- type quark sectors, respectively. In the next-to-leading order, (X15X2) scalars are 
coupled with the up-type quark sector, and Xi4 contributes to the charged lepton and down- 
type quark sectors, and then the mass ratio of the electron and down quark is reproduced 
properly. We also add Z4 symmetry to obtain relevant couplings. In order to get the natural 
hierarchy among quark and lepton masses, the Froggatt-Nielsen mechanism |99] is introduced 
as an additional U{1)fn flavor symmetry, where denotes the Froggatt-Nielsen flavon. The 
particle assignments of SU{5), S4, Z4, and U{1)pn are summarized in Table 1. The U{l)p]\f 
charges i, m, and n will be determined phenomenologically. 

We can now write down the superpotential respecting 5*4, Z4, and U{l)p]\f symmetries 
in terms of the 5*4 cutoff scale A, and the U{l)p]\f cutoff scale A. The SU{5) invariant 
superpotential of the Yukawa sector up to the linear terms of Xi = 1? ■ ■ ■ ? 13) is given as 



w = yliTuT^) ® T3 ® (xi, X2) ® i^s/A + y^T^ ®T^®H^ 

+ yUK. n;) ® {N^ n;) ® (xs, X4) ® e^^^-vA^-" + mn^^ ® k 

+ y?{Nl, ® (Fi, F2, Fs) ® (xs, Xe, Xt) ® ^^5 ® e'"/(AA™) 

+ y^N^ ® (Fi, F2, Fs) ® (X5, X6, Xt) ® H,/A 

+ yi{Fu F2, F3) ® (Ti, T2) ® (X8, X9, Xio) ® ^^45 ® e7(AA0 

+ y2{Fu F2, F3) ® T3 ® (xii, X12, X13) ® (1) 

where y'^, y2, yi, |/^, , y2, yi, and y2 are Yukawa couplings. The U{l)pN charges i, m, 
and n are integers, and satisfy the conditions m — n < and 2m — n > 0. In our numerical 
study, we fix £ = m = 1 and n = 2 phenomenologically as seen in Appendix B. Then, some 
couplings are forbidden in the superpotential. We discuss the feature of the quark and lepton 
mass matrices and flavor mixing based on this superpotential. However, we will take into 
account the next-to-leading couplings as to Xi the numerical study of the flavor mixing 
and CP violation. 

2.2 Lepton sector 

We begin to discuss the lepton sector of the superpotential w. Denoting Higgs doublets as hu 
and hd, the superpotential of the Yukawa sector respecting the S4 x Z4 x U{1)fn- symmetry 
is given for charged leptons as 



wi = -3yi 



c c 
y|(^A.X9 - IrXlo) + -^(-2/eX8 + ^ + ^rXlo) 



/i45e7(AA^ 



+ y2T%leXll + ^A.Xl2 + lrXn)hd/A- (2) 
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For right-handed Majorana neutrinos, the superpotential is given as 

+ y2 [(iVe + N;n::)xz + {KK - N;NI,)x4\ e^-VA^— + MN';N';, (3) 
and for Dirac neutrino Yukawa couphngs, the superpotential is 



yz J 

+ |/fiV;(/eX5 + /^X6 + te)/i«/A. (4) 



Higgs doublets hd and gauge singlet scalars 9, Xi ci-re assumed to develop their VEVs as 
follows: 

{K) = Vu, (hd) = Vd, (/i45) = (0) = 0, 

((X3, X4)) = {U3, U4), ((X5, X6, Xr)) = («5, Uq, U7), 

((X8, X9, Xw)) = Ug, Uio), ((Xll, X12, X13)) = {Uii, U^, U13) , (5) 

which are supposed to be real. Then, we obtain the mass matrix for charged leptons as 

/ ag/V2 -aio/V2\ / \ 

Ml = -3yi\%5 -2as/VQ ag/VQ aio/VQ + y2Vd 0, (6) 
V / ai2 "13/ 

while the right-handed Majorana neutrino mass matrix is given as 

'A2™-"(?/fA"A + |/2^a4A) |/fA2"^-"a3A 0\ 

Mn=\ yfA^'^-^agA A2"*-"(i/fA"A-?/fa4A) . (7) 

MJ 

Because of the condition m — n < 0, (1,3), (2,3), (3,1), and (3,3) elements of the right- 
handed Majorana neutrino mass matrix vanish. These are so called SUSY zeros. The Dirac 
mass matrix of neutrinos is 

^2^5/76 -ae/V6 -oj/Vq'^ 
Md = y?X^Vu I «6/V2 -«7/v^ I +y2Vu I I , (8) 


where we denote = ui/ k and A = 9/k. 

In order to get the left-handed mixing of charged leptons, we investigate M^Mi. If we 
can take vacuum alignment {us,ug,uio) = (0,^9,0) and (mh, M12, Mis) = (0,0,^13), that is 
«8 = ctio = an = «i2 = 0, we obtain 

/O -3yi\^agV45/V2 \ 
= -3yiA^«9M45/V6 , (9) 
\0 y2ai3VdJ 
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then MjMi is as follows: 



mJmi 









(10) 



where we replace 1/1^45 with yiVd- We find = — 



0, where 9l- denote left-handed 



Ql _Ql 

mixing angles to diagonalize the charged lepton mass matrix. Then, charged lepton masses 
are 







m, 



I |2 2 2 
1 2/2 I "iS^d • 



It is remarkable that the electron mass vanishes. We will discuss the electron mass as well 
as the down quark mass in the next-to-leading order. 

Taking vacuum alignment (m3,M4) = (0,M4) and {u5,Uq,U7) = {u^,u^,u^) in Eqs. ([7]) and 
([S]), the right-handed Major ana mass matrix of neutrinos turns to 



N 



'A2™-"(?/fA"A + i/2^a4A) 0_ 

A2"^-"(yfA"A-?/2^a4A) 



(12) 











and the Dirac mass matrix of neutrinos becomes 



M 



D 



'la^l^ -a^/y/Q -tts/Ve^ 
ag/v^ -a5/V2 







tts as as 



(13) 



By using the seesaw mechanism My = MJ^Mj^^Md, the left-handed Majorana neutrino 
mass matrix is written as 




3 3 3 



(14) 



where 



m\2 



m\2 



M ' A2™-"(yf A"A + ?/^a4A)' 

The neutrino mass matrix is decomposed as 



A2™-"(2/f A"A - y^a^A) 



■ (15) 



b + c 




(16) 



which gives the tri-bimaximal mixing matrix f/tri-bi and mass eigenvalues as follows: 

\ 

, nil = b , 'm-z = 3a , m^, = c . (17) 



tri-bi 



/ ^ A. 

' V6 ^/3 

]_ j_ _i_ 
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The next-to-leading terms of the superpotential are important to predict the deviation 
from the tri-bimaximal mixing of leptons, especially, Ues- The relevant superpotential in the 
charged lepton sector is given at the next-to-leading order as 

Awi = i/A,(Ti,T2) ® (Fi,F2,F3) ® (xi,X2) ® (xii, X12, X13) ® ^sM' 
+ yA,(Ti, T2) ® (Fi, F2, F3) ® (X5, X6, Xt) ® Xi4 ® i^s/A' 
+ yA.(Ti, T2) ® (Fi, F2, F3) ® (xi, X2) ® (X5, X6, X7) ® 
+ yA,{Ti, T2) ® (Fi, F2, Fa) ® (xn, X12, Xia) ® Xm ® i^45/A' 

+ yAj3 ® {Fu F2, F3) ® (X5, X6, X7) ® (X8, X9, Xio) ® ^^5 ® /A' 

+ VAfTs ® (Fi, F2, F3) ® (X8, X9, Xio) ® (Xii, X12, X13) ® ^45 ® /A' . (18) 



By using this superpotential, we obtain the charged lepton mass matrix as follows: 



Ml 



2 + ^12 ei3 

— + 622 £23 

+ 633; 



where and rrir are given in Eq. ( ITTi) . and e^j's are calculated by using Eq. ( ITSl) as 

en = Z/Ai,a5ai4Wd - 3?/Ac,«ia5Wd, 



ei2 



ei3 



1 

2' 



^(v^-i)yA.,-7(v^+i; 



yAe 



v^3 1 I 1 

— {VS- l)yAa, + t(V3 + l)yA„2 ^aiai3 - T^yA^ttsau 



^(v^+l)yA., - J(v^-l)yA. 



-2/Ad«i3ai4 



e2i 

622 

e23 



— ?/A6«5«14^'d 



7(v^-l)yA., + ^(v^+l 



ai«i3 



V3 



yAi,«5ai4 



-3 



-^(v^-l)2/A„, + ^(v^+l)yA. 

l)yAe, - ^(VS - l)yAe2| ttltts ~ ^VA^Oiisa 



Vd 



14 



Vd, 



e3i = -yA.a^agVd - SyAjOgai^Vd, 

e33 = Z/Ae«5«9^'d- 



(19) 



(20) 



Magnitudes of e^j's are of (9(q;^), where a is a linear combination of Oj's. 

Then, M/M; is given in terms of e^j, which give the non- vanishing electron mass, as 
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follows: 



M}Mi 



2\y oeii + 621 jm^ 



UV^eis + €23] 



m,, 



(21) 



Thus, the charged lepton mass matrix is not diagonal due to next-to- leading terms e^j, which 
give the non- vanishing electron mass. Since we have = 0{Xa), ttIt- = 0{a), and e^j = 
0{me), mixing angles 6*5^2' ^13 ^23 given as 



0[s = Oi^ 



^^3 = 



(22) 



Therefore, the charged lepton mixing matrix is written as 



( 1 



o 



o 
o 



( \ 



menifj 



(23) 



Now, the lepton mixing matrix U is deviated from the tri-bimaximal mixing as follows: 

U = t/^f/tri-bi- 



(24) 



The lepton mixing matrix elements Ue3, Ue2i and U^^ are given as 

'^^^"^TiO-fe))' '^"^"^TiO-K^)) ■ 

(25) 

Thus, the deviation from the tri-bimaximal mixing is lower than (9(0.01), which is rather 
small. 

Let us discuss the electron and down quark masses. The determinant of M/M/ is 



det 



MjMi 



,2_2 



- o"^/."^r ( gCii - —611621 + -621 



1 

3^ 



(26) 



where 6jj's are taken to be real for simplicity. Then the electron mass is given as 



2 



3 fl 



611621 + -621 
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2 22 



2 2 



vl- (27) 



In the same way, the down quark mass, which is discussed in subsection 2.3, is obtained as 



d 2 



(28) 
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In order to get the ratio mj : = 1 : 9, we require the following condition: 

5 {VSyA,^ - VaJ 2 {VSyA,^ - VaJ 

ai4 = ^ —ai, or = ^ —ai. (29) 

Va, Va, 

Thus, the flavon xu is introduced in our model to explain the proper ratio of the electron 
mass and the down quark mass although those masses appear at the next-to-leading order. 

Hereafter, we fix £ = 1, m = 1, and n = 2 as Frogatt-Nielsen charges, which are given in 
Appendix B. The superpotential of the next-to-leading order for Majorana neutrinos is 

Awn = yl (iVe^ ® (iV,^ Np ® (xi, X2) ® Xm/A 

+ VaM:, N') ® ® (X5, X6, Xr) ® (Xii, X12, X13) ® e/(AA) 



i2 ^ e ' /X 

+ vlsiK, ® ® (X8, X9, Xio) ® (X8, X9, Xio) ® e/(AA) 

^4 



+ yl,K ® K ® (X8, X9, xio) ® (x8, X9, xio)/a. (30) 



The dominant matrix elements of Majorana neutrinos at the next-to-leading order are written 
as follows: 

AMn = Ax 

/ ?/Ai«1«14 ?/Ai«1«14 -^Z/A2«5ai3 + ^^Aa-^ttgA 

Z/Ai"1«14 -2/Ai«1«14 -^Z/A2«5ai3 + ^I/As^Q • (31) 

\-^?/A2"5«13 + ^1/A3«9 -^1/A2«5ai3 + ^l/A3«9 1/A4«9 / 

Then, the Ues is estimated as 

1/A,«1«14 s . . 

1/2 «4 

We also consider the Dirac neutrino mass matrix. The superpotential at the next-to- 
leading order for Dirac neutrinos is given as 

AwD = y^iK, AT^) ® (Fi, F2, Fs) ® (X8, X9, Xio) ® (xii, X12, X13) ® ^^5 ® e/(A2A) . (33) 

The dominant matrix elements of the Dirac neutrinos at the next-to-leading order are written 
as follows: 

AMd = I y^XagarsVu * * \ ■ (34) 



Then, we can estimate Ues as 



Ue3 C'(a) . (35) 



Thus, the contribution of the next-to-leading terms on 11^3 is of 0{a) in the neutrino sector 
while that is 0{me/m^^) in the charged lepton sector. Therefore, it is concluded that the 
deviation from the tri-bimaximal mixing mainly comes from the neutrino sector. 
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Figure 1: The allowed region on 
sin^ 6'i2-sin plane. 
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Figure 2: The allowed region on 
sin^ 26'23-sin 6'i3 plane. 



Let us discuss the deviation from the tri-bimaximal mixing numerically. In order to obtain 
numerical result, all Yukawa couplings at the leading order are complex, and these absolute 
values and phases are taken to be randomly chosen from to 1 and —it to it, respectively. 
Other Yukawa couplings at the next-to-leading order are to be real and randomly chosen 
from —1 to 1 since the contribution from these phases on the CP violation is very small. 
Parameters a^'s are fixed as seen in section 3. We also take M = 10^^ GeV. The parameter 
ai4 is constrained to reproduce the proper ratio of the electron mass and down quark mass 
as seen in Eqs. flTTI) and f l25]) . Yukawa couplings are also constrained to give the absolute 
value of the electron mass. 

We present numerical result of the deviation from the tri-bimaximal mixing with scat- 
tering plots. Here, we neglect the renomalization effect of the neutrino mass matrix be- 
cause we suppose the normal hierarchy of neutrino masses and take small tan/3 (= 3). 



0< 300 
o 

S 200 

3 

Z 100 



We show our prediction of sin ^13 versus sin^ 9^ 
in Figure 1, where ^j/s are lepton mixing angles in 
the usual convention. That is, sin ^13 = jf/esl- In 
Figure 2, we show the prediction of sin 6*13 on the 
sin^ 26'23-sin^ 6*13 plane. The predicted upper bound 
of sin ^13 could be larger than 0.1 by tuning Yukawa 
couplings. In the case that sin ^13 is larger than 
0.15, the value of sin^ 9i2 considerably deviates from 
the tri-maximum value 1/3, that is sin^ ^12 > 0.34. 
However, the predicted points of sin ^13 distribute 
mainly in the region lower than 0.07. The mixing 
angle 623 is mainly predicted near the maximal mix- Figure 
ing angle 7r/4. The predicted points larger than 0.99 sin 6*13. 
for sin^ 2^23 cover 99%. 

We investigate precisely the predicted value 
sin6'i3 in our model. Let us calculate the expectation value of ^13. In our calculation, 
one million parameter sets are generated randomly. The number of allowed parameter sets 
in experimental constraints is 1442. These points have been plotted in Figures 1 and 2. In 
Figure 3, we present the distribution of the plot versus sin6'i3. By using this result, we can 




3: Number of points versus 
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calculate the mean value of 613 and the standard deviation. The mean value of sin 613 is 0.023 
and the standard deviation is 0.028. Thus, the expected value of ^13 is small as expected of 
0{a) in Eqs. §^ and §5^. 

It is noted that the main contribution on sin ^13 
comes from the next-to- leading term in Major ana 
mass matrix AM^r in Eq. fl31l) . which is a few times 
larger than the correction AM/) of Eq. flM|) . 

We can estimate the leptonic CP violating mea- 
sure Jcp since Yukawa couplings are taken to be 
complex. In Figure 4, we show the Jcp versus 
sin ^13. The upper bound of Jcp could be larger 
than 0.01, which may encourage the measurement 
of the CP violation in the future neutrino oscilla- 
tion experiments. However, the probable value of 
Jcp is much smaller than 0.01. We have estimated 
the mean value and its standard deviation of Jcp- 
The predicted mean value of Jcp is 2.1 x 10"'^ and 
its standard deviation is 4.5 x 10"'^. 



0.00 
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0.10 
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Figure 4: The allowed region on sin ^13- 
Jcp plane. 



2.3 Quark sector 

Let us discuss the quark sector. For down-type quarks, we can write the superpotential as 
follows: 



1 1 
— (s^X9 - &'Xio)gi + -^{-2(1^X8 + s^X9 + b^Xw)(l2 



Wd = yi 



Since the vacuum alignment is fixed in the lepton sector as seen in Eq. 
quark mass matrix at the leading order is given as 



/i45e7(AA0 

(36) 

I, the down-type 





Md = Vd\ yi\'a^lV2 yiX^a^/^Q 

2/2ai3> 



where we denote yiVd = yiv^^. Then, we have 










This matrix can be diagonalized by the orthogonal matrix U^'^ as 



(37) 



(38) 



U. 



(0) 



cos 60° sin 60° 0^ 
-sin 60° cos 60° 
1, 



(39) 
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The down-type quark masses are given as 

ml = , m1 = ^lyiA^agl^w^ , ml = \y2\'^al2Va , (40) 

which correspond to ones of charged lepton masses in Eq. (fTTI) . The down quark mass vanishes 
as well as the electron mass, however tiny masses appear at the next-to-leading order. 
The down-type quark mass matrix including the next-to-leading order is 

^11 621 ^31 

4^ + ei2 ^ + 622 632 | , (41) 

ei3 £23 mb + e-ss^ 

where q/s are given by replacing y/\. with —y/^J?) {i = ci,C2,d,f) in Eq. fl20l) . and and 
rrif, are given in Eq. fHU]) . 

In order to get the left-handed mixing, we estimate mJm^ as 



^2 ' ^ 11 -^J^ ~ 13 

i^m, + ei2)(im, + 6*22) + ene^i + £136^3 \^ + ^22? + \^2i\^ + \^23? £23^6 | ' (^2) 

By rotating the matrix M^Md with the mixing matrix uf^ in Eq. (EH]), we have 

/ C(mrfm,) i(^^3 - v^e;3)mA 

t/f MtM,f/f ^ 0(m,m.) \{V?>e\, + 6l3)m, . (43) 

\|(ei3 - y3e23)"^b 1(^/3613 €23)^6 / 

Then, we get mixing angles 6*^2, 6*13, d^z ^'^^ mass matrix of Eq. (l43l) as 

- O f = C (0.05) , ^f, = O f = O (0.005) , ei, = o('^ =0 (0.005) ,(44) 

where CP violating phases are neglected. 

Let us discuss the up-type quark sector. The superpotential respecting S^x Z4X U{1)fn 
is given as 

Wu = yi [{u^Xi + c^X2)g3 + t%qiXi + 92X2)] K/A + y^f'qsK ■ (45) 
We denote their VEVs as follows: 

((Xi,X2)) = (ni,n2) . (46) 

Then, we obtain the mass matrix for up-type quarks is given as 

/ y'^ai 
Mu = Vui y'ia2 \ . (47) 
\y^ai y'la2 1/2 

11 



The next-to-leading terms of the superpotential are also important for the prediction of 
the CP violation in the quark sector. The relevant superpotential is given at the next-to- 
leading order as 

+ yljs ® T3 ® (X8, X9, Xio) ® (X8, X9, Xio) ® H^/A^- (48) 
Then, the next-to-leading mass matrix becomes 

2/A„,aia2 yl^S'^l + (4)-ylJal~al) + ylal^ , (49) 

ylal) 

which is added to the leading up- type quark one of Eq. (147|) . In order to get the realistic 
quark mixing, we take the alignment 

«! = ^2, (50) 
then, the mass matrix of the up-type quarks is 

Mu = vA yl^al 22/X^^«2 + |/X,«?4 y>i • (51) 

After rotating by the orthogonal matrix L^^"'' as 

/ cos 45° sin 45° 0\ 
f/W = -sin 45° cos 45° , (52) 

V 01; 

we get 

({'^yla^-yla^Wi + y\oili o o 

Mu = uiMjj^ = Vu{ "0 ' {'^yl^^+ylJ<^l + ylA ^y>i 

\ V2y\a^ yl^y\al' 

(53) 

We take a phase convention in which (1, 1) and (3, 3) elements in Eq. (153|) are real. It 
is found that the magnitude of (2, 2) element is much smaller than that of (2, 3) and (3, 3) 
elements. In the limit of neglecting the (2, 2) element, the mass matrix is taken to be 
real since other phases can be removed by the phase matrix P 



(54) 
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The matrix is diagonalized by the orthogonal transformation as Vj'MuVu, where 
K = I u n I , rc = xl : , n 




rric + mt 

in which mass eigenvalues of up-type quarks are given as 



rric + mt 



(55) 



2/2 



rrir 



2/2 



m ~ y2Vu. 



(56) 



Now we can discuss the CKM matrix. Mixing matrices of up- and down-type quarks are 
summarized as 



cos 45° sin 45° 
-sin 45° cos 45° 
1 



cos 60° sin 60° 0^ 
f/rf ~ [ - sin 60° cos 60° 
1, 




ad 
-^12 - 

od I nd nd 
-<^13 + ^12^^23 





















n 










nd 


I- 


^^3 


?) 





(57) 



Therefore, the CKM matrix at the GUT scale can be written as 



10 
V' = UlUd ~ I Ti -r. 



,0 




cos 15° sin 15° 0^ 
-sin 15° cos 15° 
1, 



nd nd nd 

'^\2 '^13^^23 

ad I nd nd 
-^13 + ^12^^23 



ad 
^12 



1 ^23 
ad ad ad i 
-^23 - ^12^13 J- 



(5^ 



In the following, we suppose that 9^21 9 13, and are real. Then, nine CKM matrix elements 
at the GUT scale are expressed as 



cos 15° 



^2 



^f3^2'3)sinl5°, 



^ e% COS 15° + sin 15°, 



^12 ' 



V;° ^^f3Cosl5° + ^2'3sinl5°, 





—rte^^ sin 15° — 


0'i2 + ^^3^23)6^" cos 15° + r,(ef3 - e',20l3) 


cs 


-ri^^2e'''sinl5° + 


ne'Pcosl5° + r,{ei, + e'^20i3) , 




-rt^f3e^''sinl5° + 


^^^36*'' cos 15° - re , 


Vt'd-^ 


—Tc sin 15°e*'' — Vc 


{Ou + ^^3^23)6*" cos 15° + ^(-^^3 + e'l2e'2s 




-r^ef2sml5°e'P + 


r,e^''cosl5°-n(^2'3 + ^f2^23) , 




-rc9f2sml5°e'P + 


r^e^^e'P cos 15° + n . 



(59) 
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Putting typical masses at the GUT scale tUu = 1.04 x 10^^ GeV, rric = 302 x 10"^ GeV, 
and rrit = 129 GeV |100] . we can write CKM matrix elements in terms of p and 6fj at the GUT 
scale. We should take account the renormalization effect in order to get the CKM matrix 
elements at the electroweak (EW) scale. We use a simple formula of the renormalization for 
the CKM matrix in Ref. |101] . Then, the CKM matrix at the GUT scale [101] becomes 



v^, K° v^jhit)] , (60) 

at the EW scale. In the case of the GUT scale 10^^ GeV, we have h{t) ^ 1.05. Putting 
central values of the observed CKM matrix elements in the particle data group |103j . i.e. 
\Vus\ = 0.2257, \Vub\ = 0.00359, \Vcb\ = 0.0415, and \Vtd\ = 0.00874, we obtain a parameter 

set 

p = 123°, = -0.0340, = 0.00626, 9^^ = -0.00880, (61) 

which reproduce the experimental data. These magnitudes of 6',^'s are consistent with the 
ones in Eq. (jSj). 

In terms of a phase p, we can also estimate the magnitude of the CP violation. Let us 
calculate the CP violation measure, Jarlskog invariant Jcp |102j . which is given as 

\Jcp\ = |Im {VusV:,VubV:,} I ~ 3.06 x 10"^ , (62) 

where p = 123° is taken. Our prediction is consistent with the experimental value Jcp = 

Next, we calculate three angles of the unitarity triangle, a(or 02), /3(or 0i), and 7(or ^3), 



Putting p = 123°, we obtain a = 89.4°, sin 2/3 = 0.693 (/3 = 21.9°), and 7 = 68.7°, which 
are consistent with experimental values a = (88^5)°, sin 2/3 = 0.681 ± 0.025, and 7 = 

(77ti)° m^. 



3 Magnitudes of VEVs and Alignment 

As seen in previous section, we need relevant vacuum alignment to get the tri-bimaximal 
mixing of leptons and Cabibbo angle of quarks. The alignment of VEVs is summarized as 

(Xi,X2) = (l,l), (X3,X4) = (0,1), 

iX5, Xe, Xr) = (1, 1, 1), (X8, X9, Xio) = (0, 1, 0), (xii, X12, Xis) = (0, 0, 1), (64) 

where these magnitudes are given in arbitrary units. 

Magnitudes of = (xj)/A are determined when the quark and lepton masses are put, 
except for ai4, which appears at the next-to- leading order. Here, we have fixed £ = 1, m = 1, 
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and n = 2 as Frogatt-Nielsen charges, which are discussed in Appendix B. Then, these are 
given as 





2 


yAao 


Vu 



as = as = aio = an = oiu = 0, 
ai = a2 — 



^4 

ag 



{yfX){m3 — mi)m2M 



Qy2 1)2'^ mim^.K 



0^5 = ag = aj 



m,. 



VQ\yi\>^Vd 



"13 



y2Vd 



(65) 



where masses of quarks and leptons are given at the GUT scale. 

Putting typical values of quark masses at the GUT scale |100] . M = 10^^ GeV, A = 0.1, 
and tan/3 = 3 [vd — 55 GeV, Vu — 165 GeV) with of order 1 for absolute values of Yukawa 
couplings, we have 



ai ~ 
ag ~ 



3 X 10" 
5 X 10" 



0^4 ~ 



10"^ 
' 2 X 10" 



tts ~ 10 



(66) 



Therefore, the magnitudes of all VEVs are almost of order 10""^. Hereafter, we denote the 
averaged value as a. 

We can generate the vacuum alignment through F-terms by coupling flavons to driving 
fields, which carry the R charge +2 under f/(l)_R symmetry. 





(Xl,X2) 




(X5,X6,X7) 


(X8,X9,Xlo) 


(Xll,Xl2,Xl3) 


Xl4 


SU{5) 


1 


1 


1 


1 


1 


1 


S4 


2 


2 


3' 


3 


3 


1 




—i 


1 


—i 


-1 


i 


i 


U{1)fn 


-i 


—n 











-i 


U{1)r 























(Xl5,Xl6,Xl7) 




xl 


xt 


ixlxt) 


SU{5) 


1 


1 


1 


1 


1 




3 


1 


1 


1 


2 




-1 


-1 


i 


-1 


—i 


U{1)fn 


—z 


2£ + n 





2£ 


z 


mR 





2 


2 


2 


2 



Table 2: Assignments of SU{5), S'4, Z4, U{1)fn, and U{1)r representations for fiavons and 
driving fields. 

Three 5*4 singlets x?, X2; X3 aiid one S'4 doublet (X4, xt) for driving fields are required 
to obtain relevant vacuum alignment in our model. Moreover, an S'4 triplet (xi5, Xie? Xn) is 
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introduced as additional flavon^. Assignments of flavons and driving fields are summarized 
in Table H 

The SU{5) X 5*4 X Z4 X U{1)fn x U{1)r invariant superpotential is given as 

w' = Ki (xi, X2) ® (Xi, X2) ® (X3, X4) ® X?/A 
+ Vi iXs, X9, Xio) ® (Xii, X12, X13) ® X2 
+ V2 iXi, X2) ® (Xi, X2) ® X3 + ^3X14 ® Xl4 ® X3 

+ ^4 (X5, X6, X?) ® (Xl5, XI6, X17) ® (X45 Xs) ' (67) 

which is rewritten as 

w' = f^i [2x1X2X3 + (x? - X2) X4] X?M + ^71 (XsXii + X9X12 + X10X13) X2 



+ [V2ixl + xl) + ^3Xm] xi + -^Va (X6Xi6 - X7X17) X4 
+ -^'Hi (-2X5X15 + X6X16 + X7X17) X5 • 



Then the scalar potential is given as 



V 



X [2x1X2X3 + (Xi - X2) X4] + \rii (xsXii + X9X12 + XioXi3)r 



A 



+ h2(x^ +X2) +^3Xur + 



^74 (xeXie - X7X17) 



+ 



-^Vi (-2X5X15 + X6X16 + X7X17) 



Therefore conditions to realize the potential minimum {V = 0) are given as 

«i [2X1X2X3 + {xl - xl) X4] /A = 0, 
Vi (XsXii + X9X12 + X10X13) = 0, 



^2(Xl + X2) +^3X 



14 



0, 



(XeXie - X7X17) = 0, 
Vi (-2X5X15 + X6X16 + X7X17) = 



(6^ 



(69) 



(70) 



where Xi's are regarded as VEVs. One of the solution which satisfies these conditions is 
obtained as 



Xl = X2, X3 = 0, X5 = X6 = X7, X8 = Xio = Xii = X 



12 



0, 



2 2?72 2 

Xl4 = Xl, Xl5 = XI6 = Xl7 

V3 



(71) 



Therefore we obtain the desired alignment of VEVs in Eq. ( iMll . Next-to-leading couplings 
of flavons and driving fields could shift these alignments. Detail discussions are given in 
Appendix C. 



^As far as z 3> 1, (xi5i X161 X17) do not disturb the result in section 2. 
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4 Soft SUSY breaking terms 



We have already discussed SUSY breaking terms i.e., sfermion masses and scalar trilinear 
couplings in the flavor model [1041 1105] , the ^4 flavor model [19], and the A (54) flavor 
model [106(1107] . In this section, we study SUSY breaking terms in the framework of S4XZ4X 
U {1)fn- We consider the gravity mediation within the framework of supergravity theory. We 
assume that non-vanishing F-terms of gauge and flavor singlet (moduli) flelds Z and gauge 
singlet flelds Xi (^ = ■ ■ ■ > 14) contribute to the SUSY breaking. Their F-components are 
written as 

p'^k = _e^i^*'=^" (^djW + , (72) 

where K denotes the Kahler potential, Kjj denotes second derivatives by flelds, i.e. Kjj = 
djdjK and K^'' is its inverse. Here the flelds $/t correspond to the moduli flelds Z and gauge 
singlet flelds Xi- The VEVs of F^p^/^k are estimated as {F^^/^^) = 0(1713/2), where 171^/2 
denotes the gravitino mass, which is obtained as 7713/2 = (e^^'^^^pW/Mp) . 



4.1 Soft scalar masses of slepton sector 

First, let us study soft scalar masses. Within the framework of supergravity theory, soft 
scalar mass squared is obtained as |108] 

mjjKjj = ml/2Kij + iF^^l'd^^d^-^Ku - \F''''\'d^-^KjLd^^KMjK''^- (73) 

The invar iance under the 5*4 x Z4 x U{1)fn flavor symmetry as well as the gauge invariance 
requires the following form of the Kahler potential as 

K = Z(^)($) J2 l^il' + 4?^*^) E l^'l' + ZiS(^^\^r\\ (74) 

i=e,fj,,T i=e,^ 

at the lowest level, where Z^^\^) and Z^^^^^^{^) are arbitrary functions of the singlet flelds 
$. By use of Eq. ( |73ll with the Kahler potential in Eq. ( Tfll) . we obtain the following matrix 
form of soft scalar masses squared for left-handed and right-handed charged sleptons. 



(m|),, = I mi I , (m|),,- = | mj,^,^ | . (75) 

That is, three left-handed slepton masses are degenerate, and two right-handed slepton masses 
are degenerate. These predictions would be obvious because the left-handed sleptons form 
a triplet of 5*4, and the right-handed slepton form a doublet and a singlet of 5*4. These 
predictions hold exactly before 5*4 x Z4 x U{1)fn is broken, but its breaking gives next-to- 
leading terms in the slepton mass matrices. 

Next, we study effects due to 5*4 x Z4 x U{1)fn breaking by Xi- That is, we estimate 
corrections to the Kahler potential including Xi- Since each VEV is taken as the same order, 
the breaking scale can be characterized by the average of VEVs, such as aA. 






f ^1(1) 








,ml)ij = j 





2 







V 





2 
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In our model, the right-handed charged leptons i?p are assigned to 2 and its con- 
jugate representation is itself 2. Similarly, the left-handed charged leptons (Lg, L^, Lr) are 
assigned to 3 and its conjugation is 3. Therefore, for left-handed sector, higher dimensional 
terms are given as 

AKl = J2 ^aI mLe, L,, Lr) ® {L:, L^, L^:) ® {x^, X^+l) ® ix'i, X?+i)/A' 

1=1,3 



i=5,8,ll 

+ zgH-f )(^e, L„ Lr) ® {LI, L^, LI) ® xi4 ® Xm/A' 

+ zi";) (<i>)(Le, L^, L.) ® (l:, l^, l^) ® e ® evAl 



(76) 



For example, higher dimensional terms including (xi, X2) and (xs, Xe, Xt) are explicitly writ- 
ten as 



AK 



[X1>X5 



A2 



A2 



(77) 



When we take into account the corrections from all XiX] to the Kahler potential, the soft 
scalar masses squared for left-handed charged sleptons have the following corrections. 



■i/2) 



3/2) 



0{a'ml/^) ml + 0(aX/2) 



0{a^ml/^) 
0{d?mli^) 



(78) 



3/27 C'(«^"^3/2) "^i + C(a^m2/2) 



where a is a linear combination of aj's. 

For right-handed sector, higher dimensional terms are given as 



t.Kn = J2 zi^lmRe,R,) ® (i?:,i?P ® (x.,Xm) ® (x-,Xm)/A' 

j=l,3 

+ E R,) ® (i?:, /?^) ® (X., X.+i, X.+2) ® (X-, X^+l, X^+2)/A' 

i=5,8,ll 

+ Z2^($)(i?e, R,) ® (i?:, ® Xl4 ® Xm/A' 

+ 4?(^)(^e, R,) (^K® (Xi, X2)/A^ + zf^mRl, Rl) ®Rr® (x?, X2)/A' 
+ ($)i?. ® i?^ ® (x., X.+i) ® (X^, X?+i)/A^ 

i=l,3 

+ Zfl{mr®K®{XuX^+l,X^+2)®{xlxU^,X'i+2)/^' 
i=5,8,ll 

+ zfl{^)Rr ®Rl® xi4 ® Xm/A' 

+ z2)($)(i?e, i?^) ® {K. ® e ® evA^ 
+ zi^^ ($)i?, ® ® e ® evA^. 



(79) 
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In the same way, right-handed charged sleptons can be written as 



'm 



K) 



R(l) 



'3/2 V 
2 ^ 



"^3/2) 



3/2J 

;t.2™2 



(80) 



In order to estimate the magnitude of FCNC, we move to super-CKM basis by diagonal- 
izing the charged lepton mass matrix including next-to-leading terms. For the left-handed 
slepton mass matrix, we obtain 



2 ^^{SCKM) 
L'ij 



/ml + C(f m2 



3/2^ 



'3/2 J 



\ 



3/2) 
3/2J 



,2 I ^0/'al^2 
.A2' 



"312) 



where we take the left-handed mixing 



3/2/ 
3/2) 



nil 



3/2 



-j-a 
-a 



^.2 



-a 



12) 



In Eq. fl82!) . the unitarity is satisfied up to 0{a 



;:.2\ 



For the right-handed slepton mass matrix, we obtain 



2 \ (SCKM) 



0(aim2/2) 

where we take the right-handed mixing as 



m 



R(l) 



0(aX/2) 
C(ai"^3/2) 



3/2.' 




sin 15° 0' 
cos 15° 
1, 



X 



1 



C(«l"^3/2) 

■;^,2„ 



R{2) 



1 



In Eq. (l84l) . the unitarity is satisfied up to 0(q;^). 

Off-diagonal entries of (^l)!/*"^^^^ and 
iments |109] . Our model predicts 



^3) 



(84) 



are constrained by the FCNC exper- 



(A 



(m? 



LL 12 



2\{SCKM) 



0{a 



x2\ 



(A 



RR) 12 



(m|) 



0(a 



;?2\ 



^5) 



11 



where we take ml 



m 



as (Aii)?^P, (A 



cxp 

RR) 12 



R{1) 



m 



3/2 



. The fi ^ e'-/ experiment |110] constrains these values 



< C(10-3) [inn], when 



m-. 



m 



R 



100 GeV. On the other hand. 



the parameter space in the previous section corresponds to a ~ 10~^ and gives {All)i2, 
{Ajiji)i2 < (9(10~^). Thus, our parameter region is favored from the viewpoint of the FCNC 
constraint. 
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4.2 A-term of slepton sector 

Here, let us study scalar trilinear couplings, i.e. the so called A-terms. The A-terms among 
left-handed and right-handed sleptons and Higgs scalar fields are obtained in the gravity 
mediation as |1U8] 

hijLjRiHx = ^ hfpj^LjRjHK + h\^]^LjRjHK, (86) 



where 



K=5, 45 

hfjl^LjRjHK = -{yLjK)LjRjHKF''''K'^~'^d^,Kij (87) 
-{yiMK)LjRjH,F''^K^''^d^^KMj 
-{yijK)LjRjHKF'^^K"'^d^^KH^, 



and Khi^ denotes the Kahler metric of Hk- In addition, yuK denotes effective Yukawa 
couplings, and it corresponds to 

/ ag/V2 -aw/V2\ / 

yiJK = -3yi i-2a8/VQ ag/VQ aio/y/Q +1/2 

\ / "12 «13y 

Since the vacuum alignment indicates aio = ai2 = ctis = = 0, we get 

yiJK = -3yi I ag/VQ I + 1/2 I O I , (89) 



.0 



then 



'0 F"VV2 0^ 

hZ^K = -"-^[0 F-^/VQ 0|+^|0 1, (90) 






^0 

























^0 












where F"' = F^'/aj and F"VA = C'(^3/2)- 

Next, we estimate h\^l^. When we neglect correction terms and use the lowest level of 
Kahler potential, we obtain 

hfA = yiJKiAf + A^j), (91) 

where we estimate Af = A2 = = F"^/(ajA) ^ (9(7723/2). The magnitudes of Af = A2 
and A^ are also O {1713/2) ■ 

Furthermore, we should take into account next-to-leading terms of the Kahler poten- 
tial including Xi- These correction terms appear all entries so that their magnitudes are 
suppressed in 0{a) compared with the leading term. Then, we obtain 

(mlnh ^ ^3/2 O (a'v,) ^ O {a^v,) | . (92) 
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Moving to the super-CKM basis, we have 



\SCKM 



LRJij 



m3/2 I O (a^d) Oim^) O (a^Vd) 

'0{m,) 0{m,) 0{m,)' 
m3/2 I O (me) C(m^) O (me) 
,0(me) C(me) C(m,), 



(93) 



The FCNC measure 



^^2 



(A 



m 



LR 12 



LRJ 12 



m 



3/2 



0(me) 

^3/2 



(94) 



is predicted to be of order 5 x 10~^ for m3/2 = 100 GeV. 

This ratio is marginal compared with the bound (9(10^^) by the /i — )■ 67 experiments when 
the slepton mass is 100 GeV. Therefore, we expect that the // — 67 process will observed in 
the near future if slepton mass is 100 GeV. 



4.3 FCNC from neutrino sector 

Next, we consider the contribution on FCNC from the neutrino sector. The typical process 
of the FCNC is the /i — )■ 67 |110] . The key quantity is the magnitude of (Y^Yd) 12, which is 
given in the following matrix: 



-yf\' 



d2 



3^/2 



3^2^^ 



2 +2yfV 



3^2^^ 
3.2^^ 

3^2^' - 



- yf>^' \ 
-yfx' 

2yfxy 



yf 























«5 


«5 



(95) 



The FCNC measure on — 67 is calculated as follows |11H I112[ 1113] : 
(Am 



(A 



LL)12 



'12 
^^-'SUSY 



6ml 

167r2M|usy 



(>'r^>D)i2ln — ~ -^yfann—^ 6 x 10^^ , 



(96) 



where we put mo = Msusy, as = 10"^ A = 10^^ GeV, M = 10^^ GeV. It is concluded that 
the contribution on /i -> 67 from the neutrino sector is much smaller than the experimental 
bound iALL)lf < C(10-=^) [T09] . 



5 Summary 

We have presented a flavor model with the 6*4 symmetry to unify quarks and leptons in the 
framework of the SU (5) SUSY GUT. Three generations of 5-plets in SU (5) are assigned to 3 
of 5*4 while the first and second generations of 10-plets in SU{5) are assigned to 2 of S4, and 
the third generation of 10-plet is assigned to 1 of 6*4. These assignments of 6*4 for 5 and 10 
lead to the completely different structure of quark and lepton mass matrices. Right-handed 
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neutrinos, which are SU{5) gauge singlets, are also assigned to 2 for the first and second 
generations and 1' for the third generation. These assignments realize the tri-bimaximal 
mixing of neutrino flavors. The vacuum alignment of scalars is also required to realize the 
tri-bimaximal mixing of neutrino flavors. Our model predicts the quark mixing as well as 
the tri-bimaximal mixing of leptons. Especially, the Cabibbo angle is predicted to be around 
15°. Our model is consistent with observed CKM mixing angles and CP violation as well 
as the non- vanishing Ue^ of the neutrino flavor mixing. The deviation from 15° in \V^g\ is 
given by 0{mfi/ms). Therefore, we can adjust one parameter at the next-to-leading order to 
reproduce the observed Cabibbo angle. The non-vanishing Ue^ of the neutrino flavor mixing 
is also predicted to be ~ 0.02. 

We have also studied SUSY breaking terms. In our model, three families of left-handed 
slepton masses are degenerate and two right-handed sleptons are degenerate. Even although 
we take into account corrections due to the flavor symmetry breaking, our model leads to 
marginal values of FCNC's compared with the present experimental bounds. Therefore, we 
expect the observation of the /i — j- 67 process in the near future. 
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Appendix 



A Multiplication rule of S4 

The 5*4 group has 24 distinct elements and irreducible representations 1, 1', 2, 3, and 3'. The 
multiplication rule depends on the basis. One can see its basis dependence in our review [9] . 
In this appendix, we present the multiplication rule, which is used in this paper: 

= {aibi + a2b2)i © (-0162 + aafeOi' © l^^l^ ^ ^^^^ ) (97) 
















as) 







^aibi — 02^2, 

0261 \ / aibi 

-l{V3aib2 + a2b2) \ © |(y3a2&2 - ai&2) | (98) 
KVSaibs - 02^3) / 3 \-|(^/3a2^3 + 0163) 

aibi \ I as^i 

i(V3a262-ai&2) © -|(y3ai&2 + 02^2) | (99) 
\{\^a2bs + 0163)/ \ \{y/?>aibs - 02^3 



2 



{aibi + a2&2 + 0363)1 © 1 / )f > , > , > 

' 7g(-2ai0i + 0262 + 0363 

3 

02^3 + 0362^ 

«l63 + «3^1 I © I 0163 - 0361 I (100) 
,0162 + 02^1 







("A 


0.2 


1 «l 




\as) 


3' 






^(a2&2 - a^bs) 




( ^(^2^2 — Osb 

[aibi + 0262 + 03^3)1 ® 1 / o L , L 



^(-2ai6i + 0262 + a3&, 



3 



^0263 + 03^2^ 

aibs + a-A I © I ai&3 - | (101) 
^162 + 026 




, , , /7s(2ai6i - a2&2 - a3&3 

(ai6i + a2&2 + a363)i'©(^^ _ ^^^^^ 




02^3 + ^3^2 

aibs + asbi | (102) 
0162 + 02^1, 



More details are shown in the review [9]. 
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B Determination of ^, m, and n 

The charged lepton masses in Eq. ( fTT!) give 

vQ\yi\Xh!d 

Therefore, if we take \yi\ = 1, X = 0.1, and = 6.86 x 10"^ GeV, we get 

ag ~ 5.1 X 10^-^ . (104) 

Suppose that magnitudes of all are same order 10~^ as seen in Eq. ( I66l) . So we take i = 1, 
which gives ag ~ 5.1 x 10"'^. 

We take the right-handed neutrino mass M as follows: 

M = 0(10^2) _ ^105) 



As seen in Eqs. f[T5|) and (1171) . parameters a and c should be comparable. Therefore, we have 

2/2^A^"a4A ~ M = 0(10^2) GeV . (106) 
Again suppose that magnitudes of all are same order 10^^, we get following equation as 

y^A ~ C(10^^-") . (107) 
where A = 0.1. On the other hand, m and n satisfy the condition: 

<m<n<2m . (108) 
Therefore, we have the smallest value of n which satisfies Eq. fll08j) as 

m = l, n = 2. (109) 

C Next-to-leading terms of the scalar potential 

We consider the next-to-leading terms of the scalar potential. In the next-to-leading order, 
the SU{5) X 5*4 X Z4 X U{1)fn x U{l)fi invariant operators which couple to driving fields are 
given as follows; 

• Coupled with Xi- 

(Xi,X2)®(Xi,X2)®x;®07A", 

(X5, X6, Xr) ® (X5, X6, Xr) ® X? ® e^^+VA^^^", 

(xii, X12, X13) ® (xii, X12, X13) ® X? ® e2'+"/A2^+", 

(X3, X4) ® (X5, X6, Xr) ® (X5, X6, Xt) ® X? ® e27(AA2^), 
(X3, X4) ® (Xii, X12, X13) ® (Xii, X12, X13) ® X? ® 6^7 (AA^^), 

(X5, X6, xt) ® (X8, X9, xio) ® (xii, X12, X13) ® X? ® e2^+7(AA2^+"), 

(X8, X9, Xio) ® (X8, X9, Xio) ® (X8, X9, Xio) ® X? ® e''+"/(AA2^+"). (110) 
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• Coupled with X2- 

(X5, X6, Xr) ® (X5, X6, Xr) ® (xii, X12, X13) ® xll K 

(Xii, X12, X13) ® (xii, X12, X13) ® (xii, X12, X13) ® xl/J^- 

• Coupled with Xg: 

(xii, X12, X13) ® (xii, X12, X13) ® X3 ® e^VA^^ 
(X3, X4) ® (x5, X6, xt) ® (x5, X6, xt) ®xl® e2^-"/(AA2'-"), 
(X3, X4) ® (xii, X12, X13) ® (xii, X12, X13) ®xl® e2^-"/(AA2^-"), 

(X5, X6, Xt) ® (X8, X9, Xio) ® (Xii, X12, X13) ® X3 ® e2V(AA2^), 
(Xs, X9, Xio) ® (Xs, X9, Xio) ® (Xs, X9, Xio) ®xl® 6^7 (AA'^). 

• Coupled with (X45X5): 

(Xs, X9, Xio) ® (Xii, X12, X13) ® (Xi5, X16, X17) ® (X4> X5)/A, 

(X5, X6, X?) ® (X5, X6, X?) ® (Xll, X12, X13) ® (Xl5, XI6, X17) ® (X4' X5)/A^ 
(X5, X6, Xt) ® (X8, X9, Xio) ® (X8, X9, Xio) ® ® (Xl5, Xl6, Xl?) ® (X4> X5)/A^ 

(Xii, X12, X13) ® (Xii, X12, X13) ® (Xii, X12, X13) ® (Xi5, X16, X17) ® (X4' xD/A^- 
As seen in Eqs. fl67j) and flllOp . we can write the superpotential which couples to x? as 

= Ki (xi, X2) ® (Xi, X2) ® (X3, X4) ® X?/A 

+ /^2 (Xi, X2) ® (Xi, X2) ® X? ® 0"/A" + fi:3Xi4 ® Xi4 ® X? ® 0"/A" , 
which is rewritten as 

w^o = /€i [2x1X2X3 + (x? - X2) X4] X?/A + «:2A"(x? + X2)X? + «3A"Xi4X? • 
Then, we obtain 

«^i [2x1X2X3 + (X? - X2) X4] /A + fi:2A"(x? + X2) + «^3A"X?4 = . 
Inserting X14 = "^(Xi + X2) i'^ Eq. f lTT]) into this equation, we have 



nil 



(112) 



113) 

(114) 
(115) 
(116) 



V2 



K2 - K^— A" + fi:iX4/A 

^73 



Taking X3 — 0, we get 



xl 



K2 - fi;3— ) A" - K1X4/A 



X2 + 2fi;iXiX2X3/A = 



xi 

xl 



1^2 



) A^ 



:ii7) 



:ii8) 



As far as ^104 is much larger than (^2 — K32|)A", the vacuum alignment xi = X2 is guaranteed 
approximately. Therefore, magnitude of deviation of the vacuum alignment is parameter 
dependent. Calculating other next-to-leading operators, we can find easily that other vacuum 
alignment is deviated at most of order a ~ 0.01. 
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